This program examines tracts that only partially intersect with community districts.

Load relevant files

library(sf)
## Linking to GEOS 3.9.0, GDAL 3.2.1, PROJ 7.2.1
library(tidyverse)
## -- Attaching packages --------------------------------------- tidyverse 1.3.1 --
## v ggplot2 3.3.5     v purrr   0.3.4
## v tibble  3.1.4     v dplyr   1.0.7
## v tidyr   1.1.3     v stringr 1.4.0
## v readr   2.0.1     v forcats 0.5.1
## -- Conflicts ------------------------------------------ tidyverse_conflicts() --
## x dplyr::filter() masks stats::filter()
## x dplyr::lag()    masks stats::lag()
library(magrittr)
## 
## Attaching package: 'magrittr'
## The following object is masked from 'package:purrr':
## 
##     set_names
## The following object is masked from 'package:tidyr':
## 
##     extract
path_name <- "../data/"
load(paste0(path_name, "bl", ".RData"))
glimpse(bl)
## Rows: 39,903
## Columns: 4
## $ bl_id    <chr> "200910518012010", "200910514003006", "200910515002005", "200~
## $ bl_name  <chr> "Block 2010", "Block 3006", "Block 2005", "Block 1009", "Bloc~
## $ bl_area  <dbl> 38500.823, 51711.677, 104537.047, 73933.762, 95854.528, 19934~
## $ geometry <POLYGON [°]> POLYGON ((-94.64213 38.9774..., POLYGON ((-94.63966 3~
load(paste0(path_name, "cd", ".RData"))
glimpse(cd)
## Rows: 59
## Columns: 4
## $ cd_id    <dbl> 106, 108, 113, 102, 129, 116, 114, 101, 105, 103, 107, 109, 1~
## $ cd_name  <chr> "East Side", "Old Northeast", "Greater Downtown", "Blue Valle~
## $ cd_area  <dbl> 281713850, 118237219, 181920811, 216842678, 412671242, 264282~
## $ geometry <POLYGON [°]> POLYGON ((-94.52337 39.0941..., POLYGON ((-94.50777 3~
load(paste0(path_name, "icd", ".RData"))
glimpse(icd)
## Rows: 15,748
## Columns: 19
## $ bl_id         <chr> "290950023002014", "290950063002017", "290950022003001",~
## $ bl_name       <chr> "Block 2014", "Block 2017", "Block 3001", "Block 3004", ~
## $ bl_area       <dbl> 14983.197, 9597.627, 18151.722, 18380.541, 17691.230, 17~
## $ st_id         <chr> "29", "29", "29", "29", "29", "29", "29", "29", "29", "2~
## $ co_id         <chr> "29095", "29095", "29095", "29095", "29095", "29095", "2~
## $ tr_id         <chr> "29095002300", "29095006300", "29095002200", "2909500220~
## $ bg_id         <chr> "290950023002", "290950063002", "290950022003", "2909500~
## $ people        <chr> "62", "0", "22", "26", "42", "49", "43", "0", "12", "20"~
## $ housing_units <chr> "18", "0", "9", "18", "22", "22", "14", "0", "1", "10", ~
## $ land_area     <chr> "14996", "9606", "18167", "18396", "17706", "17473", "32~
## $ water_area    <chr> "0", "0", "0", "0", "0", "0", "0", "0", "0", "0", "0", "~
## $ block_lat     <chr> "+39.0894021", "+39.0439974", "+39.0932160", "+39.093298~
## $ block_lon     <chr> "-094.5318307", "-094.5610343", "-094.5240738", "-094.52~
## $ cd_id         <dbl> 106, 106, 106, 106, 106, 106, 106, 106, 106, 106, 106, 1~
## $ cd_name       <chr> "East Side", "East Side", "East Side", "East Side", "Eas~
## $ cd_area       <dbl> 281713850, 281713850, 281713850, 281713850, 281713850, 2~
## $ geometry      <GEOMETRY [°]> POLYGON ((-94.5324 39.0887,..., POLYGON ((-94.5~
## $ in_area       <dbl> 14983.19694, 9597.62653, 18151.72239, 18380.54058, 17691~
## $ pct_in        <dbl> 100.0000000, 100.0000000, 100.0000000, 100.0000000, 100.~
load(paste0(path_name, "red", ".RData"))
glimpse(red)
## Rows: 39,903
## Columns: 11
## $ st_id         <chr> "20", "20", "20", "20", "20", "20", "20", "20", "20", "2~
## $ co_id         <chr> "20091", "20091", "20091", "20091", "20091", "20091", "2~
## $ tr_id         <chr> "20091050000", "20091050000", "20091050000", "2009105000~
## $ bg_id         <chr> "200910500001", "200910500001", "200910500001", "2009105~
## $ bl_id         <chr> "200910500001000", "200910500001001", "200910500001002",~
## $ people        <chr> "52", "47", "51", "58", "138", "68", "35", "81", "39", "~
## $ housing_units <chr> "33", "29", "27", "26", "70", "34", "18", "28", "12", "1~
## $ land_area     <chr> "27826", "26904", "26385", "27075", "74489", "58378", "2~
## $ water_area    <chr> "0", "0", "0", "0", "0", "0", "0", "0", "0", "0", "0", "~
## $ block_lat     <chr> "+39.0427053", "+39.0427081", "+39.0426993", "+39.042717~
## $ block_lon     <chr> "-094.6217603", "-094.6228101", "-094.6238297", "-094.62~
load(paste0(path_name, "bg", ".RData"))
glimpse(bg)
## Rows: 1,684
## Columns: 4
## $ bg_id    <chr> "201039819003", "201030714004", "201030711041", "201030712051~
## $ bg_name  <chr> "Block Group 3", "Block Group 4", "Block Group 1", "Block Gro~
## $ bg_area  <dbl> 984975.5, 36602869.0, 18313082.5, 26007918.3, 11448660.8, 243~
## $ geometry <POLYGON [°]> POLYGON ((-94.94756 39.3373..., POLYGON ((-95.18756 3~

Find community ids

icd                                    %>%
  tibble                               %>%
  distinct(cd_id)                      %>%
  arrange(cd_id)                       %>%
  pull(cd_id)                           -> cd_list
for (i in cd_list) {
  icd                                  %>%
    filter(cd_id==i)                   %>%
    tibble                             %>%
    distinct(bg_id)                    %>%
    arrange(bg_id)                     %>%
    pull(bg_id)                         -> bg_list
  print(
    paste(
      i, 
      paste(
        bg_list, 
        collapse=", "), 
      sep=": "))
  cd                                   %>%
    filter(cd_id==i)                    -> cd_i
  bg                                   %>%
    filter(bg_id %in% bg_list)          -> bg_i
  for (j in bg_list) {
    bg_i %>% filter(bg_id==j) -> bg_j
    ggplot(bg_i, aes())                   +
      geom_sf(
        fill=NA, 
        color="lightgreen")                +
      geom_sf(
        data=cd_i,
        aes(),
        fill="white",
        color="darkred")                   +
      geom_sf(
        data=bg_j, 
        aes(),
        fill="NA",
        color="darkgreen")                 -> map1
    plot(map1)  
  }
}
## [1] "101: 290950078022, 290950105001, 290950105002, 290950105003, 290950105004, 290950105005, 290950106001, 290950106002, 290950106003, 290950106004, 290950107021, 290950107022, 290950107023, 290950119003, 290950119004, 290950120001, 290950120002, 290950120003, 290950124003, 290950124004, 290950125011, 290950125012, 290950125021, 290950126001, 290950126002, 290950128032, 290950128043, 290950129033, 290950156001, 290950163001, 290950163002, 290950163003, 290950170001, 290950170002, 290950171004, 290950172001, 290950172002, 290950172003, 290950172004, 290950174001, 290950174002, 290950174003, 290950175001, 290950175002, 290959801011, 290959808021"

## [1] "102: 290950019001, 290950019002, 290950019003, 290950020001, 290950020002, 290950021001, 290950021002, 290950021003, 290950021004, 290950021005, 290950022001, 290950022002, 290950022003, 290950023001, 290950034001, 290950034002, 290950058011, 290950058013, 290950110013, 290950155001, 290950155002, 290950155003, 290950156001, 290950156002, 290950156003, 290950163001, 290950163003, 290950164001, 290950170001, 290950171001, 290950171002, 290950171004, 290959808021"

## [1] "103: 290950072001, 290950073011, 290950073012, 290950074001, 290950074002, 290950074003, 290950074004, 290950075002, 290950075003, 290950081002, 290950082001, 290950082002, 290950082003, 290950082004, 290950083001, 290950083002, 290950086001, 290950086002, 290950086003, 290950086004, 290950086005, 290950086006, 290950087002, 290950090003, 290950091001, 290950091002, 290950091003"

## [1] "105: 290950061001, 290950063002, 290950075001, 290950075002, 290950075003, 290950076001, 290950076002, 290950076003, 290950077001, 290950077002, 290950078021, 290950078022, 290950079001, 290950079002, 290950079003, 290950080001, 290950080002, 290950080003, 290950081001, 290950081002, 290950087001, 290950087002, 290950088001, 290950171002, 290950171003, 290950171004, 290950175002, 290959801011"

## [1] "106: 290950022001, 290950022002, 290950022003, 290950023001, 290950023002, 290950034001, 290950034002, 290950034003, 290950034004, 290950037001, 290950038001, 290950038002, 290950052001, 290950053001, 290950053002, 290950054001, 290950055001, 290950056011, 290950056012, 290950056021, 290950056022, 290950057001, 290950057002, 290950058011, 290950058012, 290950058013, 290950060001, 290950061001, 290950061002, 290950061003, 290950063001, 290950063002, 290950075002, 290950077001, 290950154013, 290950158001, 290950160001, 290950160002, 290950161001, 290950161002, 290950162001, 290950162002, 290950164001, 290950165001, 290950165002, 290950166001, 290950166002, 290950169001, 290950169002, 290950169003, 290950171001, 290950171002, 290950171003"

## [1] "107: 202090452004, 290950046001, 290950046002, 290950046003, 290950051001, 290950051002, 290950052001, 290950053002, 290950063001, 290950065001, 290950065002, 290950066001, 290950066002, 290950067001, 290950067002, 290950069001, 290950069002, 290950071001, 290950071002, 290950071003, 290950071004, 290950072001, 290950072002, 290950073011, 290950073012, 290950073021, 290950073022, 290950074001, 290950074002, 290950075002, 290950167001, 290950167002, 290950167003, 290950167004, 290950168011, 290950168012, 290950168013, 290950168021, 290950168022, 290950169001, 290950178001, 290950178002, 290950178003, 290950178004"

## [1] "108: 290950003001, 290950003003, 290950006001, 290950006002, 290950006003, 290950006004, 290950007001, 290950007002, 290950007003, 290950008001, 290950008002, 290950008003, 290950008004, 290950009001, 290950009002, 290950009003, 290950009004, 290950010001, 290950010002, 290950010003, 290950018001, 290950018002, 290950018003, 290950018004, 290950019001, 290950019002, 290950019003, 290950019004, 290950023001, 290950154011, 290950154012, 290950154013, 290950154021, 290950154022, 290950155001, 290950160001, 290950161002"

## [1] "109: 290950087001, 290950087002, 290950088001, 290950088002, 290950088003, 290950088004, 290950089001, 290950089002, 290950090001, 290950090002, 290950090003, 290950090004, 290950090005, 290950094002, 290950095001, 290950095002, 290950095003, 290950095004, 290950096001, 290950096002, 290950096003, 290959801011"

## [1] "110: 290950078022, 290950079001, 290950088004, 290950089001, 290950096002, 290950096003, 290950163002, 290950163003, 290950172002, 290950172003, 290950174003, 290950175001, 290950175002, 290959801011"

## [1] "111: 200910509004, 200910515004, 290950090003, 290950090004, 290950091001, 290950091002, 290950091003, 290950091004, 290950092001, 290950092002, 290950092003, 290950093001, 290950093002, 290950093003, 290950094001, 290950094002, 290950094003, 290950094004, 290950095002, 290950095003, 290950095004, 290950097002, 290950098001, 290950098002, 290950099001"

## [1] "112: 200910508001, 200910509001, 200910509002, 200910509004, 290950072001, 290950072002, 290950073011, 290950073021, 290950073022, 290950074003, 290950074004, 290950083001, 290950084001, 290950084002, 290950084003, 290950085001, 290950085002, 290950085003, 290950085004, 290950086003, 290950086004, 290950086005, 290950086006, 290950091003, 290950092001, 290950092003"

## [1] "113: 202090430004, 202099805001, 290470221006, 290950003001, 290950003002, 290950003003, 290950011001, 290950011002, 290950043001, 290950043002, 290950044001, 290950046001, 290950046002, 290950052001, 290950152001, 290950152002, 290950152003, 290950153001, 290950153002, 290950153003, 290950154011, 290950154012, 290950155001, 290950157011, 290950157012, 290950157021, 290950157022, 290950157023, 290950158001, 290950158002, 290950159001, 290950159002, 290950161001, 290950161002, 290950162001, 290950167003, 290950178001, 290950178004"

## [1] "114: 202099800001, 290470202012, 290470202021, 290470202022, 290470221002, 290470221006, 290950003002, 290950152002, 290950152003, 290950155001"

## [1] "116: 290470206022, 290470206024, 290470221006, 290470222003, 290950003002, 290950003003, 290950006001, 290950007002, 290950008003, 290950009001, 290950010003, 290950110011, 290950110012, 290950110013, 290950155001, 290950155002, 290950156002, 290950156003"

## [1] "117: 290470222003, 290470223013, 290470223023, 290950151001, 290950155002"

## [1] "118: 290470205001, 290470205002, 290470205003, 290470205004, 290470205005, 290470205006, 290470205007, 290470206021, 290470206022, 290470206023, 290470206024, 290470206025, 290470206031, 290470206032, 290470206033, 290470206034, 290470206041, 290470206042, 290470206043, 290470206044, 290470206045, 290470221004, 290470222001, 290470222002, 290470222003, 290470223023, 290950155001"

## [1] "119: 290470202024, 290470203001, 290470203002, 290470203003, 290470203004, 290470204001, 290470204002, 290470204003, 290470205001, 290470205007, 290470206033, 290470206034, 290470206044, 290470209011, 290470209012, 290470209013, 290470209014, 290470209015, 290470209021, 290470209022, 290470210014, 290470210032, 290470210033, 290470210044, 290470211044, 290470212103, 290470213072"

## [1] "120: 290470202011, 290470202012, 290470202013, 290470202014, 290470202015, 290470202021, 290470202022, 290470202023, 290470202024, 290470203004, 290470203005, 290470204002, 290470204003, 290470205001, 290470205005, 290470221001, 290470221002, 290470221004, 290470221006"

## [1] "121: 290470202022, 290470203005, 290470210011, 290470210012, 290470211011, 290470211015, 290470211016, 290470212041, 290470212042, 290470212043, 290470212044, 290470212081, 290470212082, 290470212083, 290470212084, 290470212085, 290470212086, 290470212141, 291650300022, 291650300023, 291650300031, 291650300041, 291650300042, 291650300043, 291650300044, 291650301011, 291650301012, 291650301021, 291650301022, 291650301023, 291650301031, 291650301033, 291650301034, 291650302012, 291650302013, 291650302071, 291650302072, 291650302073, 291650302074, 291650302075, 291650302076, 291650302081, 291650302082, 291650302101, 291650302111, 291650302112, 291650302113, 291650302121, 291650302122, 291650302123, 291650302131, 291650302132, 291650302141, 291650302151, 291650302161, 291650302162, 291650303081, 291650303082, 291650303083, 291650303084"

## [1] "122: 290470208021, 290470208022, 290470208031, 290470209011, 290470209021, 290470210033, 290470211034, 290470211044, 290470212101, 290470212102, 290470212103, 290470212104, 290470212121, 290470212123, 290470213031, 290470213032, 290470213033, 290470213071, 290470213072, 290470213073, 290470213074, 290470213075, 290470213076, 290470213091, 290470213092, 290470213093, 290470213111, 290470213112, 290470213121, 290470213122, 290470213123, 290470213143, 290470218071, 290470218072, 290470218081, 290470218082, 290470218083, 290470218102, 290470219003, 290470219004"

## [1] "123: 290470211011, 290470211012, 290470211031, 290470211032, 290470211033, 290470211034, 290470212081, 290470212083, 290470212085, 290470212091, 290470212092, 290470212093, 290470212101, 290470212102, 290470212111, 290470212112, 290470212121, 290470212122, 290470212123, 290470212131, 290470212132, 290470212141, 290470213031, 290470213032, 290470213033, 290470213131, 290470213132, 290470213133, 290470213134, 290470213141, 290470213142, 290470213143, 290470213144, 290470219003, 290470219004, 291650302013, 291650302132, 291650306012"

## [1] "124: 290470212141, 290470219003, 291650302011, 291650302012, 291650302013, 291650302014, 291650302082, 291650302101, 291650302132, 291650303062, 291650303081, 291650303083, 291650303084, 291650303101, 291650306011, 291650306012, 291650306025, 291659800001"

## [1] "125: 200910516002, 200910517002, 200910517003, 290950093002, 290950095002, 290950096002, 290950096003, 290950097001, 290950097002, 290950098001, 290950098002, 290950098003, 290950098004, 290950099001, 290950099002, 290950100011, 290950100012, 290950100021, 290950100022, 290950100023, 290950101031, 290950101032, 290950101033, 290950102011, 290950102012, 290950172003, 290950173004"

## [1] "126: 290950101032, 290950102011, 290950102012, 290950102031, 290950172003, 290950174003"

## [1] "127: 290950102011, 290950102031, 290950129031, 290950129032, 290950129033, 290950129034, 290950129041, 290950129042, 290950129043, 290950129062, 290950129063, 290950129065, 290950130031, 290950130032, 290950130033, 290950130034, 290950131003, 290950144001, 290950174001, 290950174002, 290950174003"

## [1] "128: 290950128022, 290950128023, 290950129042, 290950129065, 290950136061, 290950143002, 290950144001, 290950144002, 290950176004"

## [1] "129: 290950124002, 290950124003, 290950125021, 290950125023, 290950127021, 290950127031, 290950127032, 290950128022, 290950142031, 290950142033, 290950142051, 290950142062, 290950143001, 290950143002, 290950143003, 290950143004, 290950144002, 290950145011, 290950145012, 290950145013, 290950145041, 290950145042, 290950176001, 290950176002, 290950176003, 290950176004"

## [1] "130: 200910517003, 200910533011, 200910533013, 200910533023, 290950100012, 290950100022, 290950101031, 290950101032, 290950101033, 290950101051, 290950101052, 290950101053, 290950101054, 290950102033, 290950102044, 290950134011, 290950134012, 290950134052, 290950134071, 290950134072, 290950134171, 290950134172, 290950134173, 290950134181, 290950134182, 290950134183, 290950173001, 290950173002, 290950173003, 290950173004"

## [1] "131: 290950101032, 290950102011, 290950102031, 290950102032, 290950102033, 290950102034, 290950102041, 290950102042, 290950102043, 290950102044, 290950130032, 290950131003, 290950133014, 290950134011, 290950134183, 290950174003"

## [1] "132: 290950102031, 290950102042, 290950129061, 290950129063, 290950129064, 290950129065, 290950130031, 290950130032, 290950130033, 290950131001, 290950131002, 290950131003, 290950132031, 290950132032, 290950132033, 290950132034, 290950132081, 290950132082, 290950132101, 290950132102, 290950132103, 290950133014, 290950133071, 290950133072, 290950144001, 290959883001"

## [1] "133: 290370600044, 290370603072, 290370603073, 290370603091, 290370603092, 290950129065, 290950133131, 290950134161, 290950135043, 290950136061, 290950136124, 290950136141, 290950144001, 290950181021, 290950182001, 290959883001"

## [1] "134: 200910533024, 200910534092, 200910534093, 200910534291, 200910534311, 200910534312, 290370600011, 290370600013, 290370600041, 290370600042, 290370601001, 290370601002, 290370601003, 290950134012, 290950134051, 290950134052, 290950134072, 290950134101, 290950134103, 290950134161"

## [1] "201: 202090402001, 202090409001, 202090411001, 202090419002, 202090429005, 202099800001, 202099805001, 202099809001, 290470221006, 291650301012"

## [1] "202: 202099800001, 202099805001, 202099809001, 290470221006, 290950152003"

## [1] "203: 202090401002, 202090401003, 202090402001, 202090406003, 202090407001, 202090407002, 202090409001, 202090411001, 202090412001, 202090412002, 202090413002, 202090414001, 202090414002, 202090414003, 202090416002, 202090429001, 202090429004, 202090429005, 202090445002, 202099800001"

## [1] "204: 202090411001, 202090419001, 202090419002, 202090420011, 202090429001, 202090429003, 202090429004, 202090429005, 202099800001"

## [1] "205: 202090401001, 202090401002, 202090401003, 202090405001, 202090405002, 202090406001, 202090406002, 202090406003, 202090407001, 202090412001, 202090413001, 202090413002, 202090445002"

## [1] "206: 202090416002, 202090416003, 202090419001, 202090419002, 202090420011, 202090420012, 202090420021, 202090420022, 202090421001, 202090421002, 202090421003, 202090422001, 202090422002, 202090423001, 202090423002, 202090423003, 202090424001, 202090424002, 202090424003, 202090429001, 202090429002, 202090429003, 202090429005, 202099809001, 202099812001"

## [1] "207: 202090406001, 202090406003, 202090412001, 202090413001, 202090413002, 202090413003, 202090413004, 202090415001, 202090415002, 202090416001, 202090416002, 202090416003, 202090416004, 202090422001, 202090422002, 202090423003, 202090429002, 202090439032, 202090439052, 202099812001"

## [1] "208: 202090415002, 202090419002, 202090424001, 202090424003, 202090426001, 202090426002, 202090426003, 202090430001, 202090439052, 202099805001, 202099809001, 202099812001, 202099815001"

## [1] "209: 200910500001, 200910500002, 200910500003, 200910501001, 200910501002, 202090430002, 202090430003, 202090430004, 202090433011, 202090433012, 202090433013, 202090451001, 202090451002, 202090451003, 202090452001, 202090452002, 202090452003, 202090452004, 202090452005, 290950046002, 290950071001, 290950071002, 290950153002, 290950168021, 290950168022"

## [1] "210: 202090428001, 202090428002, 202090430001, 202090430002, 202090430003, 202090430004, 202090433011, 202090433012, 202090451002, 202090452004, 202099805001, 202099809001, 202099812001, 202099815001, 290950152003, 290950153002"

## [1] "211: 200910501001, 200910502004, 200910503013, 200910522011, 202090427001, 202090427002, 202090428001, 202090428002, 202090428003, 202090430001, 202090430002, 202090433011, 202090434001, 202090434002, 202090435001, 202090435002, 202090436004, 202090437001, 202090437002, 202099815001"

## [1] "212: 202090427001, 202090428001, 202090437001, 202090438022, 202090438031, 202090439052, 202090440024, 202099812001, 202099815001"

## [1] "213: 200910522011, 200910522012, 200910522013, 200910522021, 202090436001, 202090436002, 202090436003, 202090436004, 202090436005, 202090437001, 202090437002, 202090438022, 202090438023, 202090438031, 202090438032, 202090438033, 202099815001"

## [1] "214: 200910522022, 200910523051, 200910523061, 200910523062, 202090436001, 202090438021, 202090438022, 202090438023, 202090438031, 202090440023, 202090440024, 202090447023, 202099815001"

## [1] "215: 202090439051, 202090440011, 202090440012, 202090440013, 202090440021, 202090440022, 202090440023, 202090440024, 202090447021, 202090447023, 202090447042"

## [1] "216: 202090438022, 202090439052, 202090440013, 202090440021, 202090440022, 202090440023, 202090440024, 202090447023"

## [1] "217: 202090405002, 202090413001, 202090413004, 202090415001, 202090415002, 202090439031, 202090439032, 202090439041, 202090439042, 202090439051, 202090439052, 202090440013, 202090440024, 202090441011, 202090441041, 202090444003, 202099812001, 202099815001"

## [1] "218: 202090439041, 202090439051, 202090440011, 202090440012, 202090440013, 202090441011, 202090441021, 202090441022, 202090441031, 202090441032, 202090441033, 202090441041, 202090441042, 202090441043, 202090442013, 202090442021, 202090447021, 202090447031, 202090447042"

## [1] "219: 202090401003, 202090402001, 202090405001, 202090443013, 202090443021, 202090444001, 202090444002, 202090445001, 202090445002, 202090446012, 202090446031"

## [1] "220: 202090442021, 202090442022, 202090443011, 202090443012, 202090443013, 202090443022, 202090446011, 202090446012, 202090446014, 202090446021, 202090446031"

## [1] "221: 202090402001, 202090445001, 202090445002, 202090446011, 202090446012, 202090446021, 202090446031, 202099800001, 291650301012, 291650301032, 291650303064"

## [1] "222: 202090405001, 202090405002, 202090406001, 202090413001, 202090439031, 202090439032, 202090439042, 202090443021, 202090443031, 202090444001, 202090444002, 202090444003, 202090445002"

## [1] "223: 202090439041, 202090439042, 202090442022, 202090443011, 202090443012, 202090443013, 202090443021, 202090443022, 202090443031, 202090443032, 202090446012"

## [1] "224: 202090439041, 202090441031, 202090441032, 202090441041, 202090441042, 202090442011, 202090442012, 202090442013, 202090442014, 202090442021, 202090442022, 202090442023, 202090443022, 202090443032, 202090447041, 202090447042"

## [1] "225: 202090442011, 202090442012, 202090442021, 202090446011, 202090446013, 202090446014, 202090446021, 202090447031, 202090447032, 202090447041, 202090447042, 202090448041, 202090448042, 202090448043"

## [1] "226: 201030711023, 201030711041, 201030712021, 201030712042, 201030712051, 202090446021, 202090447021, 202090447031, 202090447032, 202090447033, 202090448031, 202090448032, 202090448033, 202090448034, 202090448041, 202090448042, 202090448043"

## [1] "227: 201030711023, 202090446011, 202090446014, 202090446021, 202090446031, 202090448041, 202090448043, 291650303064, 291650304014"